Quantum Process Tomography on vibrational states of atoms in an Optical Lattice 
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Quantum process tomography is used to fully characterize the evolution of the quantum vibra- 
tional state of atoms. Rubidium atoms are trapped in a shallow optical lattice supporting only two 
vibrational states, which we charcterize by reconstructing the 2x2 density matrix. Repeating this 
process for a complete set of inputs allows us to completely characterize both the system's intrinsic 
decoherence and resonant coupling. 
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In recent years, there have been remarkable advances 
in directly controlling and observing the dynamics of in- 
dividual quantum systems in a variety of domains. This 
degree of control of microscopic systems is one of the 
technological advances underlying the myriad propos- 
als for realistic quantum- information processing systems. 
For instance, a number of quantum computing propos- 
als rely on atoms trapped in an optical lattice 0, 111 > 8' 
system in vifhich a great deal of work has investigated 
coherent centre-of-mass motion^, 0, 1^, full characteri- 
sation of spin states jQ], loading of individual atoms into 
lattice sites B , , and coherent interactions between 
atoms [lOl 1 1 il l 1 2*1 . Here we demonstrate a technique for 
completely reconstructing the quantum state of motion 
of an atom trapped in a lattice well. By performing 
this density-matrix reconstruction for a complete set of 
input states, we are able to completely characterize the 
quantum evolution of the system (the "superoperator" ) , 
including decoherence. 

The development of a quantum computer will rely on 
the success of error-correction to reduce errors to an ac- 
ceptable levelfT^ 113. ITsl H^. A quantum computer is 
extremely vulnerable to errors and decoherence. Ex- 
perimental characterization of both the decoherence and 
the operations will be required in order to implement 
quantum error correction 0, 0, 20] . An arbitrary 
operation may be characterized using quantum process 
tomography 113 > or QPT. The resuh of QPT is the 
superoperator, a positive, linear map from density ma- 
trices to density matrices, which governs the evolution of 
the density matrix for the operation. Unlike a propa- 
gator, the superoperator allows for non-unitary evolution 
of the system, thoroughly characterizing decoherence, re- 
laxation and loss in a system. From the superoperator 
one can determine the types of errors which occur and 
develop procedures to reduce or eliminate errors, with- 
out requiring a priori assumptions about the underlying 
physical mechanism causing the errors 23] . 

QPT has recently been demonstrated using spins in a 
NMR system ^24], the polarization of single photons ^2^ 
and a singlet state filter for photon pairs 26]. QPT is 
performed by preparing a complete set of input density 
matrices, subjecting each to the operation being tested, 
and measuring the resultant output density matrices. 
Due to the linearity of quantum mechanics, QPT of a 



process on an N dimensional system can be achieved 
by sending in iV^ linearly independent density matri- 
ces (alternatively, it has recently been shown that one 
can use a single state in a larger Hilbert space as the 
input |2l,|27|].) 

In this experiment we perform quantum process to- 
mography using the motional states of atoms trapped in 
the potential wells of a 1-D optical lattice. We examine 
processes which are independent of the electronic state 
of the atom and are only dependent upon the motional 
states of the atom. The measurements are insensitive to 
the long-range degrees of freedom and effectively trace 
over the quasimomentum in the Bloch state picture (or 
equivalently the well index in the Wannier state picture). 
We use a shallow 1-D lattice which only supports 2 bound 
bands, which we label as ground (|0)) and excited (|1)). 
The lattice is vertically oriented, causing all atoms in 
higher energy, classically unbound states to quickly fall 
out of the lattice and become spatially separated from 
atoms which remain in bound states of the lattice (the 
Landau-Zener tunneling rates from the 3 lowest energy 
bands are 3T0~^, 14.5 and 1150 per second in increasing 
order). A typical sequence, from state preparation in 
the lattice to measurement, lasts 20 ms. 

We begin by cooling and trapping Rubidium-85 atoms 
in a standard vapour cell MOT to a temperature of 7 iiK 
with an rms radius of approximately 1 mm. During the 
optical molasses stage we turn on an optical lattice in 
a vertical orientation. The optical lattice is created by 
interfering two laser beams which are detuned 30 GHz be- 
low the Rb D2 resonance at 780.03 nm. Each beam trav- 
els through an acousto-optic modulator, each of which 
is driven by a function generator, providing control of 
the relative phase between the lattice beams. By modu- 
lating the phase of one of the lasers the lattice may be 
displaced by up to a lattice spacing within one microsec- 
ond. The beams are superposed after the acousto-optic 
modulators, and co-propagate to the vacuum chamber 
with orthogonal polarizations, to reduce phase fluctua- 
tions between the beams which would impress noise on 
the lattice. The beams are separated on a polarizing 
beam splitter near the MOT. The polarization of one 
beam is then rotated such that the beams have paral- 
lel polarizations in the MOT. The beams have an an- 
gle of 50 degrees between them, creating an optical lat- 
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tice with a lattice constant of L = 0.93 microns. The 
depth of the lattice is controlled by the intensity and de- 
tuning of the beams, and is chosen to be 18 Er (where 
Er = h^/8L'^m — /i-690 Hz is the effective recoil energy of 
the lattice) at which depth it contains two bound states. 
The energy separation between the states is ?i • 27r • 5.0 
kHz. The scattering rate from the lattice beams is on 
the order of 4 Hz, which is insignificant on the timescale 
of the experiment. 

The population in each band of the lattice is 
determined by adiabatically decreasing the lattice 
potential |2a. |29| . In order to satisfy the adiabatic crite- 
rion we must decrease the potential slower than h -lO^/s^, 
whereas the fastest decrease we use is /i • 4.1 • 10^/ with 
non-adiabatic effects appearing if the turn-off is faster 
than /i • 1.4 • 10^/s^. As the depth of the wells decreases 
the energy bands gradually get closer to the top of the 
potential. Once an energy band becomes classically un- 
bound, then all the atoms in that band accelerate down- 
wards due to gravity. Since each state becomes unbound 
at a different time, each band becomes mapped into a dif- 
ferent location in space. The spatial distribution can be 
recorded on a CCD camera by fluorescence imaging. Fig- 
ure la shows a sample spatial distribution after ramping 
down the potential over a time of 45 milliseconds. Al- 
ternatively, the lattice depth may be quickly, but still 
adiabatically, lowered to a depth of 9 i?r . At this depth 
the ground state has a Landau-Zener-limited lifetime of 
250 milliseconds while the excited state has a lifetime of 
0.5 milliseconds. The excited state atoms quickly escape 
the lattice while ground state atoms remain trapped. Af- 
ter holding the depth constant for some time (typically 
about 20 ms) , the relative populations can be determined 
by fluorescence imaging, as shown in Figure lb. The 
beams creating the lattice have a Gaussian shape, with 
rms width of 3 mm, causing the lattice to be shallower 
as we move farther away from the center. The spatial 
variation of the lattice depth causes atoms at the edge 
of the lattice to tunnel out sooner than atoms near the 
center, resulting in the curved clouds seen in Figure la. 
To reduce broadening effects due to inhomogeneous well 
depths we integrate the flourescence signal only over the 
central 600 fim of the cloud where the potential varies 
slowly. 

A sample of ground state atoms is prepared by filter- 
ing out the excited state atoms from the lattice. This 
is accomplished by reducing the well depth to a depth of 
approximately 9 recoil energies, at which point only the 
ground state is bound. The potential depth is held at 
this value for 3 ms, long enough for almost all the excited 
state atoms to escape. The well depth is then adiabat- 
ically increased back to the original depth, preparing a 
sample of atoms with up to 95% occupation of the ground 
state. 

To prepare a variety of initial states, we make use of 
our ability to displace the lattice, and of the atoms' free 
evolution. Displacement of the lattice is equivalent to a 
spatial translation of the atom cloud in the lattice's ref- 



erence frame, constituting a coherent coupling between 
the energy eigenstates. In addition to coupling the two 
bound states, this induces some transitions to unbound 
states, which are lost from the lattice. Spatial transla- 
tions change the coefficients of the states as described by 
the following equations. 



|0)=^coo|0>+cio|l)+/o.s5 
|1) =^ -cla |0) +C11 |1) -f-^oss 



(1) 



The coefficients coo,cio and cn are determined by dis- 
placing the lattice with different initial state populations. 
A typical displacement used during the measurement is 
116 nm. For this displacement we measure coo,cio and 
cii to be 0.86(2), 0.50(2) and 0.53(10) respectively, close 
to the theoretical values of 0.87,0.45 and 0.63. Dur- 
ing a time period At of free evolution in the lattice, the 
ground and excited states acquire a relative phase shift 
oiiuAt, where w is the oscillation frequency in the lattice. 
Using a combination of displacement and time delay we 
can prepare superposition states with arbitrary relative 
phase. 

State tomographv|30| is performed by projecting the 
unknown state onto a set of known states. We use a set of 
non-orthogonal states {<I>i, $4} = {|0) , |1) , \9x) , \0y)} 
whereHH 16*^) = cosS |0) + sin6i |1) and lOy) = cos 6*10) - 
I sin 6* |1) . We project onto states of the form cos 9 |0) -I- 
sinO |1) by spatially displacing the trapping potential be- 
fore separating the resulting energy eigenstates. We 
choose our displacements to be small in order to have 
negligible coupling to the higher energy, unbound states. 
We use a displacement of L/8=116 nm, generating states 
with an experimental value of 9 fa 0.5 radians. The 
state \9x) can be changed into state \9y) by a time de- 
lay of a quarter period after displacement. State to- 
mography is performed by measuring the projections 
rrii = ($i|p|$i) for all states |$i). The resulting mea- 
surements {mi, ...,7714} are used to reconstruct the den- 
sity matrix. 



P = 



mi 

(m3 — i*m4)— mi cos^ O—m^ sin^ 
2 sin 6 cos 9 



(m3+i*m4) — mi cos f? — m2 sin 
2 sin 9 cos 6 

171,2 



(2) 

Process tomography is first performed on the free evo- 
lution of a quantum state in the lattice for one period. 
This allows us to completely characterize decoherence in- 
trinsic to the lattice. To perform process tomography we 
prepare a complete set of density matrices as input states. 
The four linearly independent density matrices we used 
correspond to a ground state, prepared by filtering out 
the excited state population; a coherent state with real 
coherence prepared by displacing a ground state; a coher- 
ent state with imaginary coherence prepared by adding 
a quarter-period time delay after displacement; and a 
mixed state. The mixed state can be prepared by either 
skipping the filtering step, or by preparing a coherent 
state and waiting 3 ms for it to decohere (see discussion 
below) . 
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QPT proceeds as follows: one of the input density ma- 
trices is prepared, and characterized with state tomogra- 
phy as outlined above; the same state is again prepared 
and allowed to freely evolve in the lattice for 200/is (one 
oscillation period), and then state tomography is per- 
formed on the resulting state. Figure 2 shows the pro- 
jection of each input density matrix onto the projection 
states, before and after the 'operation'. 

The super-operator £ resulting from QPT can be ex- 
pressed in a number of ways. One common form is the 
operator sum representation. 

Pout = £ (Pin) = ^ A^pinA\ (3) 

i 

where are operational elements, often called Kraus 
operators^32i, subject to the constraint ^A\Ai — I . 

The Choi matrixjsj |^ provides a straightforward pro- 
cedure to obtain experimental Kraus operators. The 
Choi matrix is defined as 

C^Y.\'){J\®£i\^){J\) (4) 

where \i) {j\ is an outer product of basis states and 
£{\i) is the super-operator acting on the matrix given 
by the outer product \i) {j\. Then 

£ (p) = Y^C^jp.j 

where Ci J — £{\i) is the i,jth 2x2 submatrix of C and 
Pij is the i,jth element of the density matrix. The eigen- 
values and eigenvectors of C_ can then be used to deter- 
mine the canonical Kraus operators, given by Ai = ^/Klki 

where Ki is the ith eigenvalue and ki is the corresponding 
eigenvector written in matrix form. The matrix \i) {j\ is 
not necessarily a density matrix, but can be written as 
a linear combination of the measured input density ma- 
trices. Using a maximum-likelihood technique, we find 
the Choi matrix which best predicts the measured out- 
put states given the measured input states. This search 
is limited to physical, i.e., completely positive, Choi ma- 
trices. The resulting Choi matrix is shown in Figure 3. 
From the Choi matrix we find that the populations are 
preserved to within experimental uncertainty while the 
coherences decay by 36 percent. 

The same data may be visually displayed using a Bloch 
sphere representation l25l l35l | . which has the advantage of 
showing how any state on the surface of the Bloch sphere 
evolves into a new state. Figure 4a shows the initial, 
undisturbed Bloch sphere before evolution in the lattice, 
and Figure 4b shows the Bloch sphere after one oscilla- 
tion. The sphere becomes prolate, by contracting toward 
the z-axis by 36%, as expected. The Kraus operators arc 
determined from the Choi matrix. The most significant 
Kraus operators are found to be Ai — 0.90/ -I- i?i and 
A2 = —0.41(72 -t- i?2 where / is the identity matrix, 



is the z Pauli matrix and i?i , R2 are small remainders 
with magnitudes bound by Tr[R\Ri] < 0.03. The other 
two Kraus operators are insignificant on the scale of our 
experimental resolution, also satisfying a similar bound. 
Kraus operators of the form / and are consistent with 
pure dephasing, as expected for from either inter-well 
tunneling or inhomogeneous broadening. 

The coherence time is thus found to be 555 ps (2.78 
periods); this is shorter than the coherence time of 2 ms 
expected based on the width of the excited band (the 
inter-well tunneling rate) and on the variation in well 
depth across the finite Gaussian lattice beams. The 
number of oscillations is however consistent with oscil- 
lations observed in other work_4., 36, j37j. It is believed 
that the observed decoherence is caused by small-scale in- 
homogeneity, such as fringes, in the lattice beams. The 
lattice beams have been spatially filtered before travel- 
ling to the MOT. Unlike earlier experiments, we use a 
two-level system, and anharmonicity is therefore not a 
factor in the observed decoherence. 

An example of an operation necessary for quantum 
information processing is a single-qubit rotation. To 
demonstrate the applicability of process tomography to 
characterizing such operations we attempt to perform 
a rotation of the Bloch sphere by a method equivalent 
to a Rabi oscillation. We sinusoidally drive the dis- 
placement of the lattice at the trap frequency, thereby 
coupling neighbouring states coherently. Process tomog- 
raphy is performed after a single period of this drive. 
The lattice displacement is driven by applying a si- 
nusoidal phase modulation with an amplitude of tt/9 
radians to one of the lattice beams. The displace- 
ment is kept small to ensure that coupling is predom- 
inantly between neighbouring states. We test both 
a sine drive, Ax{t) = XmSm^cuot), and a cosine drive, 
Ax(t) = Xm{cos{u!Qt) — 1), where Xm — 26 nm and too 
is the oscillation frequency in the lattice and the pulse 
lasts from t = Q to t — 2tt/ujo. We again find the Choi 
matrix from a maximum-likelihood model, but find the 
Bloch sphere representation to be the most intuitive. Fig- 
ures 4c and 4d show the resulting Bloch spheres, which 
have rotations of 35.5 degrees and 36.4 degrees about 
the y-axis and x-axis respectively. The Bloch spheres 
are rotated 90 degrees out of phase with one another as 
expected for driving fields 90 degrees out of phase. The 
resulting shape is close to a sphere, but the radius has 
decreased in all dimensions. In particular, the length of 
the semi-minor axis' for the sinusoidal drive is 0.69 (while 
it should be noted that in the absence of the coherent 
drive, it decayed to a value of 0.64). A simulation using 
a truncated harmonic-oscillator model predicts a rotation 
of 35.0 degrees about the y-axis. 

We have introduced a new technique for determining 
the motional quantum state of atoms in an optical lattice, 
and applied it to a demonstration of quantum process to- 
mography. We have extracted the "superoperators" fully 
characterizing the action of several different operations 
on an arbitrary state of atoms in the lattice, specifically. 
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free evolution for one period and two different resonant- 
coupling protocols. In this way, we have characterized 
the intrinsic dephasing of the system over time, and the 
effectiveness of single-qubit rotations induced by resonant 
modulation of the lattice phase. We plan to extend these 
techniques to test the Markovian approximation; to char- 
acterize and optimize bang-bang methods|38l l39j | for re- 
moving inhomogeneous-broadening effects; and to study 
the well-depth-dependence of the decoherence, investi- 
gating the role of inter-well tunneling, Wannier-Stark 
transitions, and Bloch oscillations. The procedure can 
be extended to higher-dimensional Hilbert spaces, al- 
though the number of measurements required grows ex- 
ponentially with the dimensionality. Process tomogra- 
phy should prove essential for tailoring error-correction 
protocols to the observed behaviour of particular phy sical 
realisations of quantum- information svstems |40ll4lj |; it is 
reasonable to expect that such system-by-system tailor- 



ing may prove necessary if error thresholds on the order of 
10-"* or 10-5 are ever to be reached PJllJ, [Hill. More 
generally, it is the only method to permit complete char- 
acterisation of the evolution of open quantum systems, 
and should play a central role in the toolbox for con- 
trol and study of individual quantum systems, whether 
for quantum-information processing or for other applica- 
tions. 
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Figures: 

Figure 1) a) A fluorescence image of the state popu- 
lations in a lattice obtained by adiabatic decrease of the 
lattice potential is shown . A stepwise decrease of the 
potential leads to a clearer separation of the states as 
shown in b). 

Figure 2) Matrix of measured projections. Input den- 
sity matrices (reading left to right: ground state; mixed 
state; superposition with real coherence and superposi- 
tion with imaginary coherence) are shown along the top 
while the post-selected states are listed on the side. The 
table on the left shows the projections for the input states 
while the table on the right shows the corresponding pro- 
jections after one period. Note the decreased contrast 
in the and 9y projections for the coherent states. All 
populations are unchanged to within experimental error. 

Figure 3) The Choi matrix after one oscillation in the 
lattice characterizing decoherence. The left graph shows 
the real part of the Choi matrix and the right graph shows 
the imaginary. The matrix is dominated by real com- 
ponents at the corners. The diagonal corners represent 
the mapping of populations into populations. The off- 
diagonal corners, which map coherences into coherences, 
are significantly less than one, showing a loss of coher- 
ence. The dotted lines separates the 2x2 submatrices 
of the Choi matrix. 

Figure 4) Bloch sphere representation of process to- 
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mography. a) The initial Bloch sphere representing the of coherence, c) B.S. after sine drive showing rotation 
space of pure input states. b) after 1 period of free about the y-axis. d) B.S. after cosine drive, showing 
evolution the sphere contracts horizontally due to a loss rotation about the x-axis. 
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